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INTRODUCTION 

If  one  accepts  some  of  the  most  common  types  of  laws  for  rate  of  crack 
growth  and  crack  failure,  including  the  assumption  that  incremental  growth  is 
independent  of  the  history  of  loading,  one  is  in  a  position  to  address  himself 
to  some  important  questions  of  general  behavior,  including  that  of  the  effect 
on  crack  growth  and  life  of  different  types  of  loads. 

In  particular,  we  are  concerned  in  this  paper  with  the  effect  of  load 
sequence  on  crack  growth  and  failure,  in  cases  where  the  loading  is  a 
spectrum  of  loads  of  different  amplitudes.  To  take  an  example,  let  us 
suppose  we  have  a  gun  or  other  tube,  containing  a  small  internal  crack,  to 
which  a  series  of  internal  pressures,  some  high,  some  low,  are  to  be 
applied.  We  may  ask,  will  the  ultimate  crack  growth  be  less  if  the  high 
pressures  are  applied  first,  when  the  crack  is  shortest,  than  if  they  are 
applied  last,  when  the  crack  is  longer?  And,  since  crack  length  and  stress 
may  be  simultaneously  significant  in  cataclysmic,  e.g.,  Kj^  failure,  what 
sequence  of  loads  will  give  the  shortest  life?  And  is  this  serious 
condition  likely  to  occur  under  random  loading? 

Also,  are  there  circumstances  such  that  for  growth  under  many  cycles 
of  loading,  sequence  of  loading  is  unimportant?  If  so,  what  is  the 
equivalent  steady  stress  that  gives  the  same  growth  per  cycle  of  stressing 
as  does  a  spectrum  of  loads?  What  do  we  mean  by  the  term  cycle  of  stress¬ 
ing  for  an  irregular  spectrum? 

It  is  to  such  questions  that  this  paper  is  addressed.  Our  treatment  has 
been  heuristic  rather  than  mathematically  complete;  it  is  hoped  that  a  method  of 
treatment  of  such  questions  has  been  made  sufficiently  self-evident  that  the 
reader  can  use  it  for  himself,  for  his  own  particular  problems,  if  he  so  chooses. 


CONDITION  FOR  MAXIMUM  AND  MINIMUM  CRACK  GROWTH 


The  growth  Aa,  of  a  crack  of  length  "a",  per  cycle  of  nominal  loading 
stress,  S,  is  commonly  given  by 


Aa/cycle  =  C^  Km  =  C-,  Sm  am^“  Ycm,  m  >  2 

in  regions  of  K  or  "a"  where  the  C's  are  constants  and  K  is  the  toughness 

stress'^intensity  parameter.  S  as  used  in  the  above  formula  is  usually  taken 

to  be  S  -S  .  ,  where  S  is  the  maximum  and  S  .  the  minimum  stress  in  a 

max  min  max  min 

cycle  of  stressing.  If  the  loading  spectrum  is  irregular  in  shape,  we  define 
each  rise  in  stress  as  a  cycle.  Yc  is  a  parameter,  constant  or  increasing* 
with  increase  in  "a"  and  expressing  the  geometry  involved,  such  as  1.1215 
for  an  edge  crack  or  1/^cos  ira/w  for  a  central  crack  in  a  plate  of  width  w. 

We  assume  that  the  constants  C^  and  C^  may  be  different  in  different  regions 
of  crack  length  and,  correspondingly,  that  there  may  be  sudden  transitions 
in  the  rate  of  crack  growth  as  a  crack  grows  longer.  Though  we  have  some 
exceptional  evidence  of  sudden  decreases  in  rate,  this  finding  should  be 
further  investigated,  and  we  shall  assume  here  that  the  rate  does  not 
decrease  with  crack  growth.  Furthermore,  we  restrict  ourselves  to  cases 
where  the  loading  stresses  are  tensile,  thus  avoiding  the  complication  of 
complete  crack  closure  by  compression  and  ambiguous  experimental  results. 


Thus,  for  N  stresses,  S  (=  S  -S  .  ) 

max  mm 


c 


over  N  stresses, 


S 


*For  cracks  leading  to  failure;  this  restriction  eliminates  cases  of  crack 
closure  resulting  from  crack  growth  toward  a  negatively  stressed  region. 


and  the  sum  on  the  left  must  increase,  as  will  "a"  when  the  sum  on  the 
right  increases  with  additional  cycles  of  S.* 

Thus  "a"  will  be  greater  at  any  stage  of  crack  growth,  the  larger  the 
sum  on  the  right.  And  the  latter  sum  will  obviously  be  greatest  if  the 
largest  values  of  S,  occurring  in  the  entire  spectrum  of  loading,  are 
crowded  into  it  (independent  of  order,  so  long  as  they  are  the  largest). 
Therefore,  with  the  above  law  of  crack  growth  "a"  will  be  maximum  up  to  any 
number  of  cycles,  N,  if  the  largest  occurring  values  of  S  are  used  up,  in 
computing  the  sum  up  to  the  N  stresses,  on  the  right.  And  if,  on  the  other 
hand,  the  smallest  values  of  S  are  used  first,  "a"  will  be  minimum. 

The  intermediate  cases  of  crack  growth  for  n  cycles  may  be  found,  in 
order  of  severity,  by  progressively  decreasing  the  sum  of  the  n  stresses  on 
the  right.  The  probability  of  occurrence  of  any  number  n  of  selected  loads 
out  of  a  total  of  N  occurring  in  a  representativet  cyclic  graph  of  loading  is 

p  _  n!  (N-n) ! 

K  “  N! 

The  final  value  of  "a",  i.e.,  the  total  growth  corresponding  to  use  of  . 
all  the  cycles  of  S,  is  a  fixed  value,  independent  of  the  order  of  applica¬ 
tion  of  the  stresses,  S. 

*In  conformance  with  restrictions  made  up  to  this  point,  it  will  be  observed 
that  we  confine  ourselves  to  the  case  where  the  increments  are  linearly 
additive.  That  is,  we  confine  ourselves  to  laws  of  growth  for  which  the 
incremental  growth  per  load  application  is  independent  of  the  history  of 
loading  and  environment.  Later  on  it  may  be  noted  that  we  have  implicitly 
restricted  our  failure  criteria  in  the  same  way. 

^Representative"  here  means  long  enough  to  determine  frequencies  of  occur¬ 
rence  of  all  loads.  Although  the  loading  spectrum  of  a  machine  under  test 
may  be  controlled,  in  practice  one  is  interested  in  ensemble  characteristics. 


Figure  1  illustrates  by  a  simple  example  the  above  formula  for  an 
initial  crack  length  a  =  1  and  a  loading  stress  sequence  of  three  cycles, 
S  =  1,  2,  3.  (The  stresses  and  lengths  may  be  considered  to  be  measured 


km 


in  any  units  the  reader  desires,  e.g.,  S  in  ksi,  "a"  in  inches.)  It 


demonstrates  the  method  of  computation  for  "a"  as  well  as  the  Aa/a  of  the 


above  equation,  and  shows  the  equality  of  both  sides  of  the  equation. 


Our  argument  above  will  perhaps  be  made  clearer  by  considering  the 


above  growth  law  expressed  in  integral  form.  This  representation  is 


possible  because  the  steps  Aa  are  very  small  and  numerous.  We  have 


da 

m/2  vm 


Sm  dN 


Suppose,  for  simplicity  Y  is  constant. 


if  m>2 


1  -  (aQ/a) 


m/2-1 


Y  m  am/2'1(m/2  -  1) 
c  o 


Sm  dN 


if  m=2 


in  a/ao=C2 


Here  "a"  is  the  crack  length  after  any  specified  number,  N,  of 


cycles  of  stressing.*  In  either  case,  m  >  2  or  m  =  2,  the  larger  "a"  is. 


the  larger  the  left  hand  term  is,  and  this  will  be  greatest  when  the 


integral  on  the  right  is  greatest,  that  is,  when  it  contains  the  largest 


values  of  S.  Thus  "a"  is  maximum  when  the  integral  on  the  right  contains 


the  largest  values  of  S. 


*So  long  as  "a"  remains  in  a  region  where  C 2  is  constant. 
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Figure  1.  Simple  Example.  Growth  Law  Aa/cycle  of  stressing  =  S  a. 

Crack  Length  a  =  aQ  +  E  S^a  and  E  Aa/a  =  ES^  for  N  =  1, 
2,  and  3  Cycles  of  Stressing. 


Clearly  the  order  of  application  of  these  (N)  selected  largest 
stress  cycles  is  of  no  importance  since  the  integral  on  the  right  may  be 
conceived  as  an  area  made  up  of  the  sum  of  small  rectangular  areas  of 
base  dN  and  height  Sm,  the  order  of  addition  being  inconsequential. 

Abrupt  changes  in  C  upon  reaching  certain  crack  lengths,  represent¬ 
ing  transitions  to  higher  rates  of  crack  growth  as  the  crack  grows  longer, 
do  not  invalidate  the  above  conclusion.  For,  obviously,  we  wish  the  crack 
to  have  a  maximum  length  up  to  the  transition  and  thereafter  we  may  con¬ 
sider  the  transition  point  with  its  associated  crack  length,  as  the  start 
of  a  new  problem.  Pushing  the  transition  points  closer  together,  we 
approach  the  case  of  a  continuously  increasing  C. 

Analytically,  if  C  is  an  increasing  function  of  "a",  it  will  go  under 
the  summation  of  integral  on  the  left,  as  Y^  does,  but  obviously  will  not 
invalidate  the  above  conclusions. 
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In  this  respect,  it  may  be  satisfying  to  consider  the  writer's 
crack  growth  law  with  transitions.  The  crack  growth  laws  between  transi¬ 
tions  are  like  those  just  considered  with  m=2  in  Km^Sm  am^“  but  the 
constants  are  theoretical.  The  constant  in  any  region  is  double  that  in 
the  preceding  region  and  for  constant  loading  stress,  the  crack  length  at 
any  transition  is  twice  the  length  at  the  preceding  transition,  i.e., 

a  =  2n  a  .  These  crack  growth  transitions  a  ,  N  ,  on  an  "a"  vs  N  plot 
no  b  n  n  r 

fall  on  a  hyperbola 

N.fi  .  ilSi  (1  -  5  /a) 

th 

where  a  ,  N  are  values  of  a  ,  N  . 
o  o  n  n 

"ao"  is  a  length  characteristic  of  the  material  considered.  For  some 
materials  the  usual  initial  crack  length  in  common  crack  growth  specimens 
does  not  greatly  differ  from  a  .  In  this  case  N  =  0,  of  course. 

‘th-”*.  [(i)i/3  (-f3 

where  E  is  Young's  modulus,  y  is  Poisson's  ratio  and  Y  is  the  yield  strength 
under  cyclic  loading  corresponding  to  stress*  S. 

*The  above  formula  for  C^h  is  for  the  usual  case  of  very  small  scale  yield¬ 
ing  near  the  crack  tip,  thus  for  small  S/Y  (  =  Sy/Y),  or,  more  generally, 
for  isotropic  loading  Sx  =  Sy  with  essentially  isotropic  yielding. 

A  more  generally  applicable  formula  for  the  case  where  the  loading  stress  Sx 
along  the  crack,  and  the  loading  stress  Sy  perpendicular  to  the  crack,  rise 
and  fall  together  is 

(2/33 1/3  [E/(l-u2)]'4/3  Y’2/3S  2 

cth  - - £ - 

jl  -  (VT/2)  (Sy/Y)  [(1  -  Sx/Sy)  ]  j 2/3 

(This  formula  is  the  same  as  it  was  in  the  corresponding  footnote  of  the  first 
(1977)  issue  of  this  report  except  that  a  term  then  supposed  to  render  the 
formula  applicable  to  orthotropic  yielding  has  been  removed.) 
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For  simplicity,  it  is  this  hyperbola,  which  is  periodically  correct, 
that  we  consider  as  our  growth  law,  instead  of  the  actual  law  with  transitions 
Differentiating  N, 

£ 

--  =  -  _£  (s  =  Const,  for  each  increment  and  cycle  of 

a  Lth  a“  crack  growth.) 

Hence 

(21n2)  a  Aa  v-^  ? 

- —  =  C2^S 

a 

N  stresses 

2 

where  C  =  C^/S  =  Const.  Or,  in  integral  form, 


(Note  that  m  =  4,  page  4,  also  gives  an  "a"  vs  N  hyperbola  of  this  shape, 
with  stresses  to  the  fourth  power  instead  of  second  power  as  here,  a 
difference  not  noticeable  in  constant  stress  tests  in  which  growth  rate 
constants  are  determined  empirically.) 

The  form  of  these  expressions  validates  the  conclusions  already 
drawn.  The  maximum  crack  length  up  to  any  number  of  cycles  corresponds  to 
application  of  the  cycles  of  greatest  stress,  independent  of  order. 


Figure  2  shoves  the  effect  of  load  sequence  by  extending  the  simple 
example  of  Figure  1  to  all  possible  sequences  of  the  three  loading  stresses 
given  there.  Since  there  are  three  distinct  loading  stresses,  there  are 


x 

X 


* 


six  cases  in  all,  S  =  1 ,  2,  3 ;  S  =  1 ,  3 ,  2 ;  S  =  2,  1,  3;  S  =  2,  3,  1; 

S  -  3,  1,  2;  S  =  3,  2,  1  for  the  first,  second,  and  third  cycles, 

respectively.  As  noted,  equal  crack  lengths  result  from  application  of  a 
set  of  loads  independent  of  order  of  application.  Up  to  any  point,  use  of 
the  largest  loads  results  in  greatest  "a"  and  use  of  the  smallest  loads 
results  in  least  "a". 

It  is  not  difficult,  in  some  cases,  to  satisfactorily  approximate  a 

spectrum  of  numerous  loads,  which  may  be  applied  in  any  sequence,  by  a 

few  loads  such  that  the  effects  and  frequency  of  application  of 

different  load  sequences  may  be  easily  handled  and  visualized  without 

the  aid  of  a  large  computer.  Thus,  for  example,  since  d  In  a  ~  Sm  dN  with 

Yc  =  Const.,  in  the  rate  of  crack  growth  formula,  if  there  are  ^ ,  n^  ...  nQ 

numbers  of  stresses  ...  Sq ,  respectively,  we  might  form  the  products 

n,  S1?,  n„  S™ .  n  Sm,  then  omit  any  stress  levels  for  which  any  of  these 

1122  oo 

products  are  unimportant  compared  to  the  remainder,  then  lump  the  remaining 
stress  levels  into  a  few  which  are  applied  for  the  total  number  of  cycles  con¬ 
tained  in  their  components,  e.g.,  for  stress  level  p,  of  these  few 

„m  „m 

n  S  +  n  S  +  etc. 

<,m  _  xx  x+1 _ 

p  n  +  n  +  etc. 
r  x  x+1 

applied  for  the  number  of  cycles  n^  equal  to  the  sum  in  the  denominator  of 
the  right-hand  side. 

Figure  3a  is  an  illustration  of  crack  growth  where  stresses  are  applied 
repetitively,  for  an  m  =  2  crack  growth  equation. 


1. 

Stress,  S 

o  Start 

>-  Cycle  1 

k>  Cycle  2 

ro 

0> 

r*H 

o 

X 

u 

3 

2. 

Aa=S^a 

a=a  +  £Aa 

0 

S 

0 

1 

0 

(D(l)  =  l 

1  +  1  =  2 

1 

(4) (2)=8 

,  r 

- - *-2+8=10 

3 

(9) (10)=90 

10+90=100 

2 

Aa 

0 

(D(l)=l 

(9) C  2) = 1 8 

(4) (20) =80 

a 

1 

2 

20 

100 

3. 

S 

0 

2 

1 

3 

Aa 

0 

(4) ( 1 ) =4 

(1)  (5)  =  5 

(9) (10) =90 

a 

1 

5 

10 

100 

4. 

S 

0 

2 

3 

1 

Aa 

0 

(4)(1)=4 

(9)(5)=4S 

(1)  (50)=50 

a 

1 

5 

50 

100 

5. 

S 

0 

3 

1 

2 

Aa 

0 

(9)(1)=9 

(1)(10)-10 

(4) (20) =80 

a 

1 

10 

20 

100 

6. 

S 

0 

3 

2 

1 

Aa 

0 

(9) Cl) =9 

(4) (10)=40 

(1) (50) =50 

a 

1 

10 

50 

100 

- J 

S  ' 


Figure  2.  Example:  Effect  of  Sequence  of  Loading  Stress  S  on  Crack  Growth 
Aa/ cycle  of  StressingsS^a.  Initial  (S=0)  Crack  Length,  a=l. 

S=l,  2,  and  3.  Note  Independence  of  Order  of  Application,  e.g., 
a=100  for  S=l,  2,  and  3  in  all  cases  and  a=10  for  two  cycles, 

S=1  and  2,  Cases  1  and  3. 


EQUIVALENT  STEADY  STRESS 


By  the  term  equivalent  steady  stress  applied  to  a  cracked  specimen, 
we  mean  a  cyclic  stress  of  constant  amplitude  whose  application  would 
result  in  the  same  crack  growth  as  that  due  to  any  specified  portion  of 
a  loading  spectrum  of  irregular  amplitude,  if  the  number  of  cycles  of 
load  application  were  the  same  for  each  type  of  loading. 

As  so  defined,  it  may  have  different  values  for  portions  of  the 
spectrum  made  up  of  different  loads,  but  it  will  have  a  unique  value  for 
the  whole  spectrum  since  this  includes  all  the  loads. 

The  expressions  for  equivalency  are  apparent  from  the  summation  or 
integral  expressions  we  have  already  obtained.  Since  for  both  the  constant 
amplitude  steady  and  the  irregular  stress  the  crack  growth  is  to  be  the 
same,  the  left-hand  sides  of  these  expressions  must  be  the  same  whether 
produced  by  a  steady  stress  or  an  irregular  stress.  Thus  if  there  are  N 


cycles  in  the  portion  of  the  spectrum  being  considered,  and  if  the  steady 


equivalent  stress  is  S  ,  the  sum  or  integral  on  the  right  must  be  equal 

eq 


to  S™  N  . 
eq  r 


That  is. 


,m 


eq  r 


•/ 


Sm  dN  = 


,m 


N  stresses,  S. 
r  ’ 


for  the  power  laws;  m  =  2  for  the  envelope  hyperbola. 


Thus  if  there  are  n,  stresses  S, ,  n_  stresses  S_,  ...  n  stresses  S  , 

1  12  2  r  r 


such  that 


Part  A 


With 

In  a  /a  = 

1  0 

da  =  S2  a/100 
dN 

S2  nj/100  In  a2/a1 

=  S2  n2/100 

in  a3/a2 

=  S2  n  /100 
3  3 

In  a./a  = 
2  0 

In  a./a  + 
1  0 

In  a2/a1 

In  a  /a 

3  0 

=  In  a./a  + 
2  0 

in  a3/a2 

Region,  i, 

0  (Start) 

1 

2 

3 

Cycles,  m 

0 

10 

2 

1 

S. 

0 

1 

2 

3 

1 

S? 

l 

m/100 

0 

0.10 

0.08 

0.09 

zs2 

ni/100 

0 

0.10 

0.18 

0.27 

[  a./a.  1(=s2n./100) 

1  l-l  1  1 

0 

0.10 

0.08 

0.09 

Eln  a./a  -  (=£S?n./100)  o  In  a./a  =  0.10  In  a_/a=0.18  lna,/a  =  0.27 

1  11  10  Z  O  3  o 

ai  a  =  1  a  =  1.11  a  =  1.20  a  =  1.31 

o  12  3 


Part  B 


Let  Zm  =  Nt(=13);  m/Nt  = 

in  a3/ao  =  (sf  ^  ♦  S*  „2  .  S*  „3)/100 

’  (SJ  fl  *  S2  f2  *  S3  f3>  V100  '  Seq  V100 

Thus  S2  =  S2f.  ♦  S2f-  +  s2f,  =  1  (10/13)  +  4(2/13)  +  9  (1/13)  =  27/13 

eq  1  1  2  2  3  3 

S  =1.44 

_ _ _ eq _ 

Figure  3.  Example  with  m=2,  da/dN*S2a/100.  Amplitude  (Part  A)  and  Equivalent 
Steady  Stress  (Part  B)  when  ni  Cycles  of  Sj,  n2  Cycles  of  S2,  and 
n3  Cycles  of  S3  are  applied  in  succession. 
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we  have* 


_m  ..  _  „m  cm  cm 

S  N  =  n  S  +  n_S  +  .  .  .  +n  S 
eq  r  1  1  2  2  r  r 


=  (f^J  +  f2S“  +  .  .  .  +  frS“)  Nr 

whence 

sm  =  f.  s"1  +  f_S®  +  .  .  .  +  f  sm 

eq  1  1  2  2  r  r 

in  which  m  =  2  for  the  hyberbola  envelope  law  as  well  as  the  m  =  2  power 


Figure  3b  shows  the  calculation  for  S  for  the  simple  case  of  Figure  3a, 

eq 

where  three  loads  were  applied  in  succession  for  different  numbers  of  cycles. 

For  another  example,  suppose  that  the  distribution  of  frequencies  is 

Gaussian  in  N  cycles,  o  being  the  standard  deviation, 
r 


where  S  is  the  mean  value  of  the  cyclic  stress  rise  S  *  S  -  S  .  ,  so 
M  max  min 


and  thus 


5m  -  q 
eq“  & 


S  =  J  and  thus 

£/ H  ■  4  *  GTl-fc*) 

s/M -(■*)- 


•See  Appendix  for  a  cyclic  relationship  in  lieu  of  equivalent  stress. 


where  X  =  (S-Sy)/o 

Since  the  lower  limit  (least  value)  of  S  is  zero,  the  lower  limit 

of  X  is  -  (S  / a ) . 

y 

If  this  and  the  upper  limit  are  large,  we  may  assume  the  integration  is 
between  -°°  and  + 


Let  us  assume  this  is  the  case  and  that  m  =  2.  Then 


The  exponential  function  is  positive  and  symmetrical  about  X  :  0  so 

that  the  middle  term  in  the  brackets  makes  no  contribution  and  we  have 

c2  2  _2  "Gaussian  S  " 

S  =  a  +  S  eq 

eq  V  - - -2 - 

Naturally,  if  o  is  zero,  there  is’  no  scatter  and  S  =  S 

eq  y 

The  case  of  the  parabolic  distribution,  with  its  finite  cut-off,  is 
treated  for  general  values  of  m  as  part  of  the  problem  of  probability  of 
failure. 

PROBABILITY  OF  FAILURE 


In  this  section  we  consider,  in  an  illustrative  way,  the  probability 

of  failure  of  a  piece  of  material  containing  a  crack  or  notch  and  subjected 

to  the  various  succession  of  loads  of  a  given  spectrum  of  loading.  It 

is  presumed  that  the  piece  will  fail  when  the  loading  stress  reaches  a 

critical  value,  S  .  In  current  cases  of  interest  concerning  cracks, 

this  value  of  S  is  commonly  obtained  from  an  association  of  critical  tough- 
me 

ness  value,  K^,  and  crack  length.  The  latter  depends  on  the  number  of  cycles 
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of  loading  together  with  the  equivalent  steady  stress  to  which  the  piece 


has  been  subjected.  The  practical  limitations  of  this  approach  to  Smc  and 
probability  of  failure  will  be  briefly  considered  along  with  notch  failure. 

We  first  consider  the  case  of  random  loading  and  then  a  class  of 
ordered  loading.  In  the  latter  case  we  consider,  in  particular,  the  effect 
of  the  condition  for  maximum  crack  growth. 

1.  Random  Loading 

Let  us  assume  that  the  distribution  of  both  the  maximum  stresses  S  , 

m 

S  =  S  ,  and  the  stress  rises,  S  =  S  -S  .  are  parabolic*  with  the  same 
m  max  max  mm  r 

standard  deviation  a.  Let  S  be  the  mean  value  of  S  and  S  be  the  mean 

my  max  y 

value  of  the  rises,  S. 

Thus,  for  the  maximum  stresses 

•  -°/r  5  <W  * o/r 

with  f  =  0  for  S  >  cry's  and  for  S  <  ~a/5 
m  m 

so  that  the  probability  that  S  >  S  is 
r  m  me 


f  = 


3/5 

20a 


1  - 


(S  -S  ) 
m  my' 

c  2 
5a 


0/5” 

f  3/5 

.  (W 

J  JUo 

5„2 

s  -s 

■  - 

d  (S  -S  ) 
m  my' 


me  my 


-  (S  -S  )/o  1- 

mc  my 

-  J  L 


(3/5/20)  {2/5/3  -  (S  -S  )/o  I  I  1- (1/15)  (S _ -S  T/a 


me  my 


V1 


•Note  that  the  constants  in  the  expressions  used  below  are  not  the  same  as 
that  of  the  peak  region  of  the  Gaussian  distribution.  Here  o  is  the  a  of 
the  parabolic  distribution,  not  the  a  of  the  Gaussian. 


14 


■r.% 


VV'.V-W.V  vw>\1  *  ^vv- v”'1  ."v  v*  ',"v  ;vvvv/.’  y  •.*  • 


for  -a/5  s  (S  -S  )  <  a/5, 
m  mir  ’ 


with  P  =  0  for  S  £  a/J  and  for  S  <  -a/ 5. 

m  m 


Similarly  for  the  stress  differences,  S  =  S  -S  .  , 

max  min 


3/5 

20a 


1  - 


(s-y‘ 

5a2 


,  -a/5"  <  (S-S^)  <  a/5" 


with  f  =  0  for  S  >  a/5  and  for  S  <  -a/5, 


so  that  the  equivalent  stress,  S  ,  is  given  by 


eq 


w 


0/5 

»  <■ 

(  S”N  (g) 

(s-s  y 

1  -  ; 

/  "W 

So2 

d(S-S  ) 
v 


where  the  crack  growth  per  cycle  is  proportional  to  S 
Since 


m 


sm  *  (X+S  )m  with  X  =  S 

-S,  » 

V 

V 

♦a/5 

-  - 

sm  -  ~  f  (X*S  )m 
eq  20a  1  u 

,  .  *L 

c  2 
5(7 

-a/5 

•  - 

6(m+2)  [(6*T)ra+2-K-S*T)m+2]  -  t T(6+t ) m+2_ ( - 5 +T ) m+2 ] 
«3  (m+1)  (m+2)  (m+3)  (2/3) 

where  6  *  a/5  and  t  =  S  . 


Thus*  if 


m  = 

1, 

S  = 
eq 

S 

y 

„2 

2 

m  = 

2, 

S  = 
eq 

„3 

S 

y 

_3 

0 

2 

m  = 

3, 

S  = 
eq 

_4 

S 

y 

4 

+ 

3S  o 
y 

,  J2 

m  = 

4, 

S  = 
eq 

S 

y 

+ 

6  S  a 
y 

The  form  of  the  expression  for  S  for  m  =  2  is  the  same  as  it  is  for  the 


Gaussian  distribution  which  we  derived  heretofore  even  though  the  approximate 
form  of  the  peak  region  of  the  Gaussian,  which  is  also  parabolic,  has  different 
constants. 


Knowing  Se  from  this  formula  we  can  compute  the  crack  length  after  any 


number  of  cycles  of  loading  prior  to  failure.  If  m  =  2  we  may,  for  example, 
have  the  writer's  crack  growth  law,  or  his  periodically  correct  hyperbolic 
simplification  (page  7)  of  it.  In  the  latter  case.t 


.m+2 


*  (6+t)"'*  ■  EO  +  EE  by  binomial  expansion 

where  SO  =  E  1st,  3rd,  5th,  etc.,  terms 


ZE  =  E  2nd,  4th,  6th,  etc.,  terms 
whence  For  even  m  For  odd  m 


Sm  =  3[6(m+2)E0-TEE] 

eq  63(m+l)(m+2)(m+3) 


sm  _  3[6 (m+2) EE-tZO] 

eq  63(m+l) (m+2) (m+3) 


+(1)  Let  it  take  N  cycles  to  get  to  a 
o  o 


No  "  No  =  (2  ln  2/Cth) (1  ~  ao/ao5  :  p*  6) 


(2)  Let  it  take  (N  +  Nq)  cycles  to  get  to  "a"  (changing  the  meaning  of  N 
on  page  6) 


(N  ♦  N  )  -N  *  (2  ln  2/C  .  )  (1  -  art/a) 


o'  o  '  •  "th-' '  o' 

(3)  Subtract  (1)  from  (2),  above,  to  get  the  N  cycles  it  takes  to  get  from 


a  to  "a." 
o 


N  «  -  (2  ln  2/Cth)(ao/a  -  iQ/ao) 


which  gives 


a 


V  [l  -  (a0/i0)(CthN)/(2  ln2)J  (v. :  p.  7) 


in  all  of  which  C.,  «CS‘ 


a  •  aQ  /[I  -  (aQ  /aQ) (C  Sgq‘  N)/(2  In  2)] 


where  a  is  a  material  constant  (~~2  mm),  a  is  the  measure  of  the 
o  o 

initial  crack  length  (half  length  of  a  central  crack  in  a  wide  plate)  and 

S  is  the  equivalent  range  of  loading  stress  acting  perpendicularly  to  the 
eq 

crack,  derived  above  for  the  parabolic  distribution  . 

If  a  critical  stress  intensity  factor,  K^,  is  applicable  to  the 
material  used  and  is  known  for  the  temperature  where  failure  is  to  be 
evaluated, 

KIC  -  tS„c  Yo 
S„c  -  KIC/(/a  V 

where  Smc  is  the  peak  loading  stress  corresponding  to  KIC,  is  a  parameter 
which  is  constant  or  increases  for  cracks  leading  to  failure,  and  "a"  is 
a  crack  length  or  depth  which  we  find  from  our  crack  growth  formula.  Thus 

S»c  -  '“iC-V'  Vf'-'VVl'  SeqN)/2  ln  2)1 
and  our  probability  of  failure  is  given  by  our  formula  for  P,  that  is,  the 
probability  of  occurrence  of  a  stress  equal  to  or  greater  than  S  : 

>>(s„ls.cWVs720)  j Vs7M(s>c-s„u) Wi  -  (i/i«emc-snu)2/„2i 

where  S_.  is  the  mean  value  the  peak  stresses, 
mu 

For  (S  -S  )  >  o  /5,  the  probability  P  is  zero.  Since  S  is 

the  value  of  Kj^/C/a  Y£)  it  is  not  a  loading  stress  per  se,  but  is  the 

stress  required  to  cause  fracture,  a  stress  value  which  the  loading  stress 

may  or  may  not  equal.  Thus  S  may  exceed  any  loading  stresses  applied  in 

me 

our  bounded  loading  spectrum  in  which  case  there  is  zero  probability  of  the 
occurrence  of  any  loading  stress  as  large  as  S  . 


As  the  crack  length  measure  a,  increases  with  increase  in  N,  or 

Kjc  decreases  with  (ordinarily)  decrease  in  temperature,  Smc  =  K^/i/a"  may 

decrease  from  a  safe  high  value  until  it  lies  within  the  loading  range,  in 

which  case  there  will  be  a  probability  of  failure  given  by  the  above 

formula  for  P.  If  S  falls  so  that  S  -S  «  -  a/5,  P  =  1  and  failure 

me  me  mu 

is  certain. 

In  fact  K  may  fall  so  far  and  so  rapidly  with  temperature  decrease 

on  the  lower  side  of  the  toughness  transition  temperature  that  S^c  is 

practically  independent  of  changes  in  length  a.  In  this  case  Smc  may  be 

lower  than  the  operating  range  of  stresses  and  the  failure  may  be 

regarded  as  a  stress  failure  independent  of  crack  length,  associated  with 

operation  below  the  transition  temperature,  rather  than  as  a  KI(,  failure. 

Failures  insensitive  to  crack  length,  but  responsive  to  stress  and 

transition  temperature  have  been  observed  and  studied  by  this  writer. 

On  the  other  hand,  there  are  cases  where  the  crack  growth  necessary 

to  bring  Smc  down  to  the  operating  range  of  stress  is  so  great  that  only 

a  practically  insignificant  amount  of  life  would  remain  even  if  KjC 

failure  did  not  intervene.  In  practice,  through-the-section  yielding 

may  even  occur  in  the  net  section  below  the  crack,  before  a  critical 

failure  stress  is  attained,  such  that  the  resultant  deformation  may  render 

the  piece  useless.  In  this  case,  however,  the  critical  stress  S  would 

me 

not  be  found  from  the  usual  formula  for  which  presuppose  small  scale 
yielding. 

Fatigue  life  may,  of  course,  be  estimated  by  use  of  the  equivalent 
steady  stress  if  one  employs  a  fatigue  life  law  conformable  in  type  to 
that  used  in  our  analysis  of  crack  growth.  Power  law  approximations  to 
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any  exact  law  are  suitable.  Thus  if  we  have  a  stressed  notch  to  deal 


with,  we  may  calculate  a  life  to  first  macroscopic  crack  appearance  and  a 

subsequent  life  using  for  crack  depth  a  depth  which  includes  the  original 

notch  depth.  If  appreciable  crack  growth  occurs,  the  type  of  analysis 

for  probability  of  failure  with  simple  cracks  which  we  have  illustrated 

above,  would  be  applicable.  The  probability  of  failure  prior  to  cracking 

may  be  gotten  from  the  above  formula  for  P,  omitting  the  crack  growth 

equation.  S  in  this  case  might  be  determined  from  destructive  tests 
me 

on  notched  material  made  after  a  few  cycles  of  loading,  or  estimated 
from  the  fracture  stress  F  of  the  material,  if  this  is  known,  by  using 
theory  which  takes  account  of  yielding  about  the  notch  root. 

2.  Ordered  Loading 

Here  we  deal  with  a  finite  number  of  loads  of  different  amplitudes, 
as  with  a  gun  that  may  be  loaded  with  a  variety  of  charges  and  whose  life 
is  limited  by  crack  growth  if  not  terminated  by  the  application  of  a 
sufficient  load  or  corresponding  K^.  Vie  wish  to  know  what  material 
toughness  would  be  required  to  withstand  any  selected  load  of  the  loading 
spectrum,  after  the  application  in  any  sequence,  of  any  other  selected 
loads . 

For  the  most  severe  condition,  these  loads  would  be  chosen  to  give 
maximum  crack  growth  and  the  loading  stress  applied  after  their  application 


(1)  BEEUWKES,  R. ,  Jr.  "Characteristics  of  Crack  Failure,"  Surfaces  and  Inte 
faces,  Syracuse  University  Press,  Syracuse,  NY,  Vol.  II,  1968,  p.  277. 

(2)  BEEUWKES,  R. ,  Jr.  "Determination  of  Fracture  Stress  and  Effective  Crack 
Tip  Radius  from  Toughness  (Kjq)  and  Yield  Strength  (Y),"  Army  Materials 
and  Mechanics  Research  Center,  AMMRC  TR  78-44,  October  1978. 

(3)  TRACEY,  DENNIS  M. ,  and  FREESE,  COLIN  E.  "Cyclic  Plasticity  Near  a  Crack 
like  Elliptical  Flaw,"  Mechanics  of  Materials,  Vol. I,  1982,  p.  151-159. 
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would  be  chosen  to  be  the  maximum  loading  stress  in  the  entire  spectrum  of 

loads.  Although  this  final  load  may  be  included  among  those  causing  the 

maximum  crack  growth,  this  inclusion  is  unimportant  where  there  are  many 

cycles  of  loading  since  no  appreciable  growth  occurs  in  any  one  loading. 

The  crack  depth  resulting  from  application  of  any  group  of  loads  may 

be  found  as  discussed  in  preceding  sections.  If  the  final  loading  stress  in 

the  group  is  S  ,  failure  by  the  KTr  criterion  will  occur  if 
m  il 

<Sn  Yc  i  Kic- 

assuming  is  known.  Conversely  the  toughness  necessary  to  prevent  failure 
must  be  greater  than  the  product  on  the  left-hand  side  of  this  relationship. 


APPENDIX 


Cyclic  Relationship  in  Lieu  of  Equivalent  Stress 


We  had  (top  of  page  12) 


"i  s”  + 


n„S„ 
2  2 


+n  S 
r  r 


with  E  n.  *  N 
i  r 


Let  etc*  defined  by 


S1  "fl 


S®  N 
eq  r 


c® 

S2  nf  2 


Sm  N 
eq  r 


etc. 


Then 


S®  N 


S®  N  S™  N 

S®  N  -n,  -^L-E  +  „  -Sa_E  +...+„ 

2  "t2  1  "fr 


J  «  ”  11-  - 

eq  r  1  n 


whence 


i  -  1  nfi 


of  which  a  particular  case  may  be  fatigue  provided  the  n  correspond  to 

fatigue  failures  and  S  and  N  to  the  endurance  limit. 

eq  r 


i 
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